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ABSTRACT 

The equation of drift shells, traced by the guiding center motion of charged 
particles moving in a magnetic field, is discussed in terms of Euler potentials 
a and /3. Particular attention is given to fields deviating to a limited extent from 
a dipole configuration, for which it is shown that the result is related to the 
drift-shell parameter L, with an added "shell splitting function" G r A pertur- 
bation method approximately deriving Gj is described: it leads to results simi- 
lar to those found by Pennington in his perturbation derivation of drift shells. 
The use by Pennington of a divergent expansion and the 1=0 limit of the equa- 
tions obtained are also discussed. 
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EULER POTENTIALS AND GEOMAGNETIC 


DRIFT SHELLS 


INTRODUCTION 

A drift shell (sometimes termed "magnetic shell") is defined as the surface 
traced by the guiding center of a magnetically confined charged particle, as de- 
rived from the guiding center approximation of its motion. By this approxima- 
tion, a particle generally follows magnetic field lines, and therefore such lines 
will be tangential to drift shells. 

Magnetic field lines may be compactly described by the use of Euler poten- 
tials a and /3 


B = Va x V/3 (1) 

and therefore it may be surmised that drift shells are also best expressed by 
such potentials. It is the purpose of this work to approximately derive this re- 
lationship for drift shells in the geomagnetic field. For that Held, approximate 
Euler potentials may be derived by perturbation (Stern, 1967; the notation and 
results of that work will be freely used here) from the spherical harmonic ex- 
pansion of the geomagnetic scalar potential y. A similar derivation allowing for 
external sources of the geomagnetic field is also possible; its various proper- 
ties, including the form of drift shells in that case, will be discussed in a sepa- 
rate article. 
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With one exception, previous derivations of drift shells in the geomagnetic 
field have involved extensive numerical calculation of line integrals, by means 
of large digital computers. The exception is the perturbation derivation by 
R. Pennington (1967; performed in 1960 for the Argus experiment), of which 
unfortunately only a short summary (Pennington, 1961) has appeared in the 
periodical literature. The results obtained here are fully equivalent to those 
of Pennington; their derivation, however, is through the use of Euler potentials, 
which represents the more natural approach to problems of this sort. Certain 
mathematical aspects, not treated by Pennington, are also clarified. 

FIRST-ORDER EULER POTENTIALS 

Let the geomagnetic field B be given through its scalar potential y 

B = - Vy (2) 

and let y be split (assuming tilted dipole coordinates) into a dipole component 
y 0 and a sum of higher harmonics y , assumed to be of the order of ey 0 , with 
e « 1: 


y 0 - a g° (a/r) 2 cos 6 


y i 


m~ n 

^~~' (a/r) n+1 P™(0) / g™ cos m <p + h“ sin m <p\ 


i=2 m=0 


( 3 ) 

(4) 
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In the summation of y , the n = 1, m = l terms are absent, because tilted - 
dipole coordinates are used; some of the terms there are actually small enough 
to be ignored at this level and could be relegated to a second-order component 
y 2 (Stem, 1967), but this point will not be stressed. We denote by a and fi the 
first-order Euler potentials, as derived from the dipole potentials by first order 
perturbation (Stern, 1967) 

a - + dj 

p = p* + p, 

with 

a 0 = a g° (a/r) sin 2 0 (5) 


otj - a 


m- n 

EL 

n = 2 m-0 


(a/r ) n sin 2n 6 


v^(0) + e. 




cos m cp + h™ sin m cp 


} 


( 6 ) 


A, = a9 


(?) 


m~ n 

= ^ E E 

n=2 m- 0 


(a/r) 


n" 1 


sm 


2n-2 


0 t“(<9) 


{ h ; 


cos m cp - g™ sin m cp 


} 


( 8 ) 


0 


Here V”(0) and t™(6) are trigonometrical polynomials introduced by Pennington, 
available in tabulated form (Pennington, 1961; Stern, 1965, 1967), and are 
additive constants (Stem, 1967), chosen in such a manner that the analogous 
additive constants of t™(<9) all vanish. 

THE DRIFT SHELL EQUATION 

Since field lines are tangential to drift shells, the equation of any such shell 
has the form 

f (a, P) = 0 (9) 

For a given particle, the associated drift shell and therefore also the func- 
tion f (a,/3) depends on the initial conditions of the particle's motion, e.g., its 
position and momentum at some given instant, or on suitable independent func- 
tions of these. By the adiabatic theory of guiding center motion it may be shown 
that in static cases no more than two such functions determine the shell — the 
magnetic moment of the particle and its longitudinal invariant. These, in turn, 
are functions only of the particle's mass and energy and of two quantities asso- 
ciated with the field, the "mirroring field intensity" B m at which the particle is 
reflected in its motion along field lines and the integral 

I = J (1 - B/B m ) M ds (10) 

B<B 

m 
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evaluated along a field line. It is the property of a drift shell (by which, indeed, 
it is defined) that the same value of I is obtained no matter which of the field 
lines tangential to it is chosen for integration. The entire family of shells in the 
geomagnetic field may thus be characterised by an equation of the form 

f (a, /3, I, B m ) = 0 (11) 

or 


a = G(/3, I, BJ 


( 12 ) 


Since the geomagnetic field may be regarded as a perturbed dipole field, one may 
expect that its drift shells reduce to those of the dipole field as all higher har- 
monics tend to zero. In a dipole field the variable ft, which is then given by 
Eq. (7), is absent from Eq. (11) due to the field's axial symmetry, so this equa- 
tion may be written (subscript zero referring to the zeroth order, i.e., the dipole 
case) 


Inverting, one gets 


a o = G 0 (I, B m ) 

1 = 


(13) 


(14) 


The function I o , defined here for later use, is none other than the integral of 
Eq. (10), expressed in (a 0 , 6, B m ) variables for the dipole field (see Eqs. 23-25). 
If one allows the perturbation to shrink to zero, Eq. (12) reduces to (13); there- 
fore, within first-order accuracy, Eq. (12) has the form 
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(15) 


a 0 + a i = G o (I, BJ + g; (A I, B m ) 

The function Gj appearing here is of the first order, and in evaluating it one 
may use zero-order relationships. Now, of all variables entering the last equa- 
tion, the integral I is the most difficult to handle; it is therefore advantageous 
to eliminate it from Gj by use of (14), even though this means reintroducing a 
to the right-hand side (although Eq. (14) is only correct to zeroth order, such a 
substitution is permissible in a first-order correction term). In the next section 
we shall discuss the deeper significance of this elimination — namely, that it 
allows labeling a shell by two parameters, the effects of which are of different 
orders and can be considered separately. We shall also replace a and fi in GJ 
by their zero-order parts a 0 and /3 0 ; for practical applications, when we want 
the equation of a shell in ( r, <9 , cp ) variables rather than in terms of a and / 3, 
one may substitute for them from Eqs. (5) and (7). The function obtained after 
all these changes will be denoted by unprimed G 1 , and the equation becomes 

a = G o (I, BJ + G, (a 0 , <P, BJ (16) 

THE DRIFT -SHELL PARAMETER L 

The preceding derivation bears a strong relation to the parameter L, intro- 
duced by Mcllwain (1961) for labeling drift shells. There has existed some mis- 
understanding concerning this parameter, and in what follows we shall try to re- 
solve it, at the same time defining the relationship between L and the present 
calculation. 
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Much of the above misunderstanding can be traced to the fact that there 


exist two different functions which are commonly denoted by L. On the one hand, 
there is the function L(I,B m ) of the adiabatic invariants, which is used to label 
drift shells; on the other hand there is also a function of position ("L at a point"), 
denoted here by L(r, <9, cp), which is obtained, for a given position in space, by 
substituting in L(I,B m ) the values of I and B m appropriate to particles mirroring 
at that point. The purpose of this section is to show that 

(1) L (I ,B m ) is a function naturally arising when classifying drift shells in a 
dipole field or in a perturbed dipole field, and 

(2) L (r,<9,<p) is an approximation to a certain choice of the Euler potential 
a, obtained by averaging Eq. (16) over B m . 

To demonstrate the first point, consider again the labeling of drift shells: as 
was noted, such labeling may be accomplished by means of the two parameters I 
and B m . This, however, is not the only possible choice, and two independent 
functions of I and B m may serve equally well. 

In particular, if the field is that of a magnetic dipole, it appears advantageous 
to choose as one of the labeling parameters the function G o (I, B m ), retaining, say 
B m for the other one. By Eq. (13) the shape of the magnetic drift shell then de- 
pends only on one characterizing parameter, and the entire collection of shells 
reduces to a one-parameter family, rather than a two-parameter one. 
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There still remains some freedom of choice left — instead of labeling with 
G Q ( I, B m ), one may use some function of it. In particular, one may introduce the 
function 

L(I. BJ = a g?/G o (I, BJ (17) 

so that Eq. (13) for a shell in a dipole field takes the form 

L(I, BJ = ag?/a 0 (18) 

By Eq. (5), a dipole shell crosses the equatorial plane at a fixed distance, and it 
can be seen from the last equation that the constant L associated with such a 
shell has the useful intuitive property of equaling this distance, as measured in 
earth radii. It is therefore to be identified with the function L (I, B m ) introduced 
by McELwain (1961). 

Of more interest is the case in which the field is not strictly a dipole field 
but has a perturbation added to it. To the first order the equation of a drift shell 
is then given by (16), and replacing there G o by L we get 

a = a gJ/L + G x (a 0 , fi Q , BJ ( 19 ) 

In this case the equation of a shell depends on both L and B m . Since, how- 
ever, this is a perturbed version of Eq. (18), the dependence is unequal: L ap- 
pears in a zero-order term, whereas B m enters only through the first-order 
correction G t . The function G t represents the variation between shells shar- 
ing the same value of L, and it therefore appears appropriate to call it the 


10 



(first-order) shell splitting function. Much of the present work is devoted to its 


explicit derivation. 

In this connection it should also be realized that no simple expression is 
available for either G o (I , B m ) or L (I , B m ). The only way to derive these func- 
tions is through the inverse function I o (a 0 , B m ), obtained by extracting I from 
G o (I, B m ) and shown in Eq. (14); this function is explicitely defined in equations 
(23)-(25) and is by no means a simple one. Because of this difficulty, an analyti- 
cal approximation to L (I, B m ), introduced by Mcllwain (1961, 1966), is nowadays 
generally used. 

Turning now to the second definition of L, by calculating "L at a point" one 
is actually deriving a function which to first order can be approximated by 

L(r,<9, <p) = i//(a,/3,B) = a gj/[a - Gj(a,/3,B)] . (20) 

On the right-hand side, a,/ 3 and B are all to be evaluated at the point (r, 0,<p). 
Since B enters only through a first-order correction term, the dependence on it 
is weak and may be averaged out 

L(r,<?,<P) = <^(a,AB)> aver . over B = X (a, /3) . 

Since x is function of a and /3 only, it may be introduced as an Euler poten- 
tial to replace a (Stem, 1967). In fact, if the correction term G x is altogether 
ignored, we find that L(r, 9, <p) closely approximates (a gj/a), with a defined as 
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in Eq. (6). One may even improve such an approximation by an appropriate 
choice of the constants g “ appearing in that equation, but this point will not be 
discussed. 

In summary, then, the use of L (I, B m ) for labeling drift shells naturally 
enters into the present treatment through the function G o , and is further aug- 
mented by the addition of a "shell-splitting" correction term G r If the pertur- 
bation is mild — as is the case near earth — the correction term is small and it 
is a good approximation to assume that the same value of L(I, B m ) characterises 
all particles attached to a given field line. This value can then be derived as 
L (r, 9 , cp) for some arbitrary point on the line — e.g., the point at which it meets 
the earth's surface, which is the choice used by most experimenters. 

On the other hand, the use of L (r, 5,<p ) as an Euler potential for labeling 
geomagnetic field lines (as has been practiced by many workers in the field) is 
not particularly encouraged. As has been shown, it is not entirely accurate, and 
although the inaccuracy (for the internal magnetosphere) is no greater than that 
introduced by the first-order approximation of Eq. (6), it is a basic one and is 
independent of mathematical precision. Moreover, the Euler potentials of a 
field depend only on its sources and structure, and it should not be necessary 
to involve adiabatic invariants of particles trapped in it for their derivation. 
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THE FIRST-ORDER SHELL SPLITTING FUNCTION 


In order to derive an approximate expression for Gj , Eq. (16) is rewritten 

a " Gj G 0 (I,B.). (21) 

This resembles (13), except that a 0 is replaced by (a - G^. We therefore 
get for the perturbed field, in analogy with (14) 

I = ^(a-Gj.BJ 

= - G i + °( e2 >- 


No first-order errors are committed by replacing 91 o (a, B m ) / 3 a by 
91 (a A , B ) /9a„, so that one finds, to first order 

o u m u 

I(a,AB m ) - I 0 (a,BJ 
G l (a o’ < * ) > B m> = " 9I o (a 0 ,B ra )/9a 0 * 

To evaluate the denominator, one has to differentiate 


loK.BJ = 


J [l - B o (a o ,0)/B m ]* x 0 (a 0 ,6) dd 


B <B 
o m 


( 22 ) 


(23) 


with 


B 0 = |Vy 0 | 


( 24 ) 


= gj O 0 /a i?) 3 sin 6 6 (1 + 3 cos 2 6) A 
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*0 = = r 2 B o /Py 0 /38) 



- - a (a 0 /a gj) 1 

sin 0(1+3 cos 2 8) A . 

(25) 

Using 

3B o /3a 0 = 

3 B 0 / a o 

(26) 



" X 0 / a o 

(27) 


and denoting by & m ("mirroring 0") the value of 0 corresponding to B ra in a 
dipole field, one gets 


BI o /3a 0 = - a,,' 1 (“l + B o /2B m ) (1 - B 0 /B m )’* X 0 d6 


re m 

- (a 2 g°/2a 0 2 ) (2 + B Q /B m ) (1 - B 0 /B m )" vs sin 0 (1 + 3cos 2 0)*d0. 


The integral derived here is the same as the integral K 4 defined by Penning- 
ton (1961, Eq. 11). It has an integrable singularity at the limits of integration. 

EVALUATION OF THE NUMERATOR 

The numerator in Eq. (22) consists of the first-order difference between 
two integrals. The first of these is 
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I (a, 0 , BJ 


(29) 


r 9 n 

= (1 - B/BJ* X d# 

Jn-8 

s 

a tfi = const. 


where B = |Vy| is the magnitude of the perturbed field, integration is performed 
with constant (perturbed) a and 0, and where 

X = 3s/Bi9| at/3 = r 2 B/(By/B0) | r>(p . (30) 

The integrand of I vanishes at 6 = <9 n and at 6 = rr -6 (subscripts for 
"(9-north" and "(9-south"); because the perturbed field is not necessarily sym- 
metric, 9 n and 9 s will in general be unequal, and may depend on <P. 

The second integral in the numerator is I (a, B m ), already defined in 
Eq. (23). It should be noted that here its first argument is a, not a 0 as in 
Eq. (23), so that 9 m is defined by B o (a, 0 m ), not B o (a 0 , 9 m ), beii^ equal to B m . 

The two integrals represent two functions I (a, 0, BJ and I o (a, B m ) that 
differ only slightly, and it is therefore natural to try to expand I, in some way, 
around I Q ; when such an expansion is then substituted in the enumerator, its 
zero-order terms cancel, leaving the first-order difference explicitely stated. 

In what follows the procedure for such an expansion, which is somewhat tricky, 
will be derived; the actual calculation, which is merely tedious, will not be 
given. 
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Three factors have to be considered in deriving the difference between I and 
I Q . First of all, there is the difference in limits of integration: this will be dealt 
with later. Secondly, I contains the perturbed variables B and 7, whereas only 
their dipole components B q and y 0 appear in I o . Finally, the form in which B 
and 7 are available is as functions of the spatial coordinates (r, 0, cp). By means 
of Eqs. (5) and (7), this dependence is easily transformed into a dependence on 
(a 0 , /3 0 , 6 ), but since it is a and fi, not a 0 and /3 0 , that stay constant during inte- 
gration, one has then to substitute in the above dependence 


a.Q - a ~ 

00 ~ 0 ~ 0 \ 


(31) 


and expand, thus expressing the variables in terms of (a, fi, 6). Because the cal- 
culation is to be accurate only to the first order of perturbation, the substitution 
(31) is only required in zero-order terms, while in first-order terms the dif- 
ference between (a 0 , /3 0 ) and (a, /3) may be ignored. 


As an example, consider the field intensity B (given for the dipole field in 
Eq. 24) 

B (a, fi, 0) = ! V 7 I 

= BoK>*) + ( V To ' V ?i)/B 0 
= B o (a 0 , 6) + Bj (a 0 ,/3 0 , 6) (defining Bj) 

= B q (a - a, , 9) + Bj(a,/3,0) 

(32) 

— B o (a,(9) - a 1 (BB o /3a) + B t (a, /3, 9) . 
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For convenience, we introduce a special notation for the first-order component 


AB(a,/3, 9) - - a 1 (BB Q /Ba) + B 1 (a,/3,9). 

Similarly, from the dipole relation 


one get 


r " r o ( a O> 0 ) 


r (a, 0,8) = r o (a,0) - a,(3r o /3a) 
= r o + Ar . 


(33) 


(34) 


(35) 


Finally, denoting for brevity 

= r 0 (a 0 ,£) 

one finds 

~by/ne = r 0 (a, 6) + Ar (36) 

where 

Ar = - ai (Br 0 /Ba). (37) 

For reasons which will become clear later, it is advisable to leave the difference 
arising from the square root terms for separate consideration. Expansion of all 
other terms to the first order yields 


r e « 


I (a, /3, B ) = 


1 - 


B o + AB' 

B 


= const. 



2A r 
r 

o 



X 0 (a, 6) dd. (38) 
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Let us write 


I 



- (B o + AB)/B m }* X 0 d0 



(39) 

B o /BJ w X 0 d0 + AI, . 


The quantity AIj is of the first order and may therefore be evaluated, to the 
present order of approximation, between the limits 6 m and rr-6 m . Making other 
allowable approximations, and replacing (a, /3) by (a 0 , /3 0 ) during integration, 
we get 



At the corresponding point in Pennington's calculation, the square root in 
the first part of (39) was expanded by the binomial theorem: 

{1 - (B o + AB/Bjf = (1 - B o /B m )*{ 1 - AB/(B m - B o )f 

(41) 

= C 1 ~ B 0 /B b )* - (AB/2 BJ (1 - B o /B m )^. 


In addition, the difference between the limits of the two integrals in (39) was 
ignored, leading to 


J »0 m 

AB (1 - B o /B m )-* X 0 d0 + AI X + 0(e 2 ) . (42) 

TT-6L 


18 



Actually the binomial expansion is divergent when (B m - B o ) is less than AB 
and as a result, the integrand in (42) has a singularity at its limits, where the 
quantity it tries to approximate really tends to zero. The singularity is inte- 
grable, however, and it is shown in the appendix that Eq. (42) is in fact correct 
to the first order in e. 

The explicit derivation of (I - I 0 ) from here on will not be described: as 
was noted before, it is lengthy, though not too difficult. Since the first order 
quantities AB, Ar and Ar all involve a, or 7, in a linear fashion, the final ex- 
pressions for ( I - I ), and consequently also for , split up into a sum of terms, 
each of which is proportional to one of the harmonic coefficients g” or h™. The 
calculation of these terms leads to a series of integrals, first given by Pennington 
(1961), which combine to form a set of functions a™ (Pennington's notation, re- 
tained here; not related to the Euler potential a) of the mirroring angle & m . The 
functions a" {& m ) have been evaluated numerically and tabulated by Pennington 
(1961) and more extensively and accurately by Stern (1965); using them one gets 

m~ n 

o, = <i/g°) £ 2>° /ae ? r + ^ 

n= 2 m= 0 

(43) 

cos ra <p + h™ sin m cpj. . 

Combining this with Eqs. (4), (5) and (19) gives the shell's first-order equation 
as 
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1/L - (a/ r) sin 2 6 


+ (l/g°) ^ ^( a A) n sin 2n 6 [\T (9) + a“ (0 m )] { g“ cos m cp + h” sinmcp} . 

n=2 m=0 

It should be noted that the integration constants £ “ cancel out in the final 
result. This was to be expected, since in characterizing drift shells, as in Eq.(9) 
and the ones following it, the only requirement for a and /3 was that they be con- 
served along field lines, and it was not necessary that they "match" (Stem, 1967). 

THE 1=0 OMIT 

One may test Eq. (44) by considering the limit 1=0, for which the shape of 
a drift shell (which then contracts to a line) may be separately derived from first 
principles. Such a test also checks whether the ratio of expanded integrals tends 
to the correct limit when the integration ranges become very short. 

In a dipole field, shells with 1=0 describe particles confined to the equa- 
torial plane. Since (aL) then equals their (constant) distance from the dipole, one 
finds in that case the limit of the function L (I, B m ) as 

L (°- B J = (g?/B m ) 1/3 . (45) 

Substituting this into (44) and inserting there rr/2 for the mirror angle 6 m gives 
the first-order equation of the (0, B m ) shell, by the present perturbation scheme, 
as 
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B 2 

01 


(g?) 2 (a/agj) 6 


(46) 

+ 6g° ^^(a 0 /agj) n+5 - a“ ( 77 / 2)} (g™cosmcp + h^sinmcp). 

n , m 

For an alternative derivation of the shell’s equation, we note that a particle 
with 1=0 will always be located at points at which the field intensity is B m . 

This gives to first order the condition 

B 2 = B 2 = B 2 + 2 v 7o . y 7l 

(47) 

= (g ?) 2 (a/ag °) 6 F(6) - g° ^(a 0 /agj ) n+5 Q™(£) [ g™ cos m <p + h* sinm<p] 


where 

F ( 6 ) = (1+3 cos 2 9 ) sin -12 Q 

Q^(0) = 6F(0) {\™( 9 ) + £”} 

- 2 sin“ (2n+10 > 6 {2 (n + 1) cos 6 PJJ (0) - sin 9 (dP™/d6>)} . 

To obtain the shell's equation, one would now have to eliminate 9, using the addi- 
tional condition that the position of the particle on any field line corresponds to 
the minimum of B there. However, it is sufficient to assume that the angle 
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77/2 - 6 


\ = 

between the particle's location and the equatorial plane, is of the first order in 
smallness. With this assumption one may substitute 

e = tt/2 

F(0) = 1 

not only in first order terms but also in zero-order ones, since 

3F/30| w/2 = 0 

and thus 

F(0) = F ( 77 / 2) + 0(e 2 ) . 

Comparing then (48) with (47) gives 

6{£“ - a“(7r/ 2)} = (^(77/2) = + 2(dP™/d0)U /2 (48) 

from which one finds the requirement 

a™ (tt/2) = - |(dP”/d0)L/ 2 (49) 

which indeed is met by the functions a“ (Stem, 1965, last equation). 
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APPENDIX 


PROOF OF EQUATION (42) 

Let the integration region be divided at 6 = tt/2 , and consider the contribu- 
tion to the first two terms in (39) from angles 6 smaller than this value 


AI 

n 



(B o + AB/Bjf X 0 d 9 



B o /Bj«X 0 Ad. (A-l) 


Defining a constant A <9, which may be of either sign 

A 6 > = 0 - e m 

n m 

= - AB(0 n )/(3B o /30)|* 

in v/ m 


(A -2) 


one changes variables in the first part of (39) to 

\p - 6 — A 0 . 


To the first order 


B o (a, 6) - B o (a,i p) + tSd (BB q /B^) 

/ 

AB (a, /3, &) = AB(a,/3, \jj) 

X 0 (a , 6 ) = X 0 (a.'/O + M (BX o /B0) 

giving, to the same order 


23 


r *» w 

J [i - (B o + AB)/Bj x 0 d& 


f t 

^ TT/2-&& 


1 - B o + AB + A0(3B o /B0) /B m } (X 0 + M 3X 0 /3i/») d<// 


Jim 1^ 

A^< (1 - B o /B m ) H (3X O /B0) d0 + [l - B o (a,7r/2)/B m ] X 0 (a t n/2) 


O - [B o + AB + A0 (BB 0 /B</r)]/B m J- * 0 di/> 


(A-3) 


The square root appearing in the last integral may be written 


1 - B 0 /B m )* |l - [AB + A0(3B o /B0)]/[B m - B c ]}“ . (A-4) 


Unlike in the case of Eq. (41), the second factor here may be expanded, for the 
ratio contained in it is less than unity at all points at which the integrand of (A-3) 
is real. There exists no problem at the end of the integration range, for both 
numerator and denominator vanish there. One may therefore write (reinstating 
6 as integration variable) 
r e m 

AI n = (-1/2BJ [AB + A0(3B o /30)] (1 - B o /B ra )-* X Q dd 
n/2 

re m (A-5) 

+ A e I (1 - B o /B m ) M (3X o /30) de + M [l - B o (a,77/2)/B m ] X 0 (a, n/2). 

7T/2 

However, integration by parts shows that all terms involving A 6 cancel. A 
similar cancellation occurs in the remaining half of the integration range, and 
the only terms left are those in (42). 
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Figure 1 schematically illustrates the meaning of the preceding result. 
Graphs (1) and (2) give the integrands appearing in equation (A-l), and the curve 
bounding the hatched area describes their difference. Graph (3) describes the 
binomial approximation to this curve as used in Eq. (42), an approximation which 
evidently breaks down near d = & m . In spite of this breakdown, the use of the 
approximation in evaluating integrals leads to correct results, because the area 
between graph (3) and the two orthogonal axes shown in the figure equals (to the 
order of approximation) the hatched area beneath the curve. 
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